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Abstract 

It is important to understand the mechanism of high-temperature superconductivity. It is obvious that the interaction with large 
energy scale is responsible for high critical temperature T c . The Coulomb interaction is one of candidates that bring about high- 
temperature superconductivity because its characteristic energy is of the order of eV. There have been many works for the Hubbard 
model including three-band d-p model with the on-site Coulomb repulsion to investigate a possibility of high-temperature supercon¬ 
ductivity. It is, of course, not trivial whether the on-site Coulomb interaction leads to a pairing interaction between two electrons. 
We argue that high-temperature superconductivity is possible in the strongly correlated region by using the variational Monte Carlo 
method for the two-dimensional t-U-J-V model. The exchange interaction I and the nearest-neighbour attractive interaction V 
cooperate with U and will act to enhance the critical temperature. 
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1. Introduction 

A basic question in the study of superconductivity is on 
an upper limit of the critical temperature T c . It has been sug¬ 
gested that the maximum of T c is in the range from 30K to 40K 
when the superconducting transition is based on the electron- 
phonon interaction [1, 2], The prime reason of low T c for the 
electron-phonon mechanism is that the Debye frequency a>o 
is of the order of 100K and T c is reduced by more than one 
digit than eu w . In general, the electron-phonon coupling con¬ 
stant A is small. Thus there is an upper limit for T c within the 
electron-phonon mechanism. The interaction of large energy 
scale will be responsible for high temperature superconductiv¬ 
ity. The electronic origin mechanisms of superconductivity are 
worth studying in the search of high temperature superconduc¬ 
tors. The electronic interaction comes from the Coulomb inter¬ 
action between electrons. From this point of view, it is impor¬ 
tant to consider the intra-atomic Coulomb interactions that are 
of the order of eV. Then a question arises as to whether the su¬ 
perconductivity is indeed induced by the Coulomb interaction. 
This subject has been studied intensively by using electronic 
models with the on-site Coulomb interaction U including the 
two-dimensional (2D) Hubbard model [3-26], the d-p model 
[27-33] and also the ladder Hubbard model [34-38]. 

It remains matter of controversial as to whether the 2D Hub¬ 
bard model accounts for high-temperature superconductivity. 
Most of the results of quantum Monte Carlo methods have failed 
to show the existence of superconductivity [13,14,18], and some 
results, however, still keep open the possibility of supercon¬ 
ductivity [15, 24]. We must note that quantum Monte Carlo 
methods are methods being applicable only in the weakly cor¬ 
related region where UIt is not large and is at most from 4 to 


5. The existence of superconductivity in strongly correlated re¬ 
gions should be considered further. This may account for prop¬ 
erties of high-temperature cuprate superconductors. 

2. Superconductivity in Strongly Correlated Region 

Superconductivity in the strongly correlated region is de¬ 
scribed by using the 2D Hubbard model. The Hamiltonian is 
given by 

11 = X +U Yj n ^ ,lii ■ ( 1 } 

ij i 

{ tjj } are transfer integrals and U is the on-site Coulomb energy. 
The transfer integral tij is non-zero f, ; = t for nearest-neighbor 
pair (ij) and tij = t for next-nearest neighbor pair (ij); other¬ 
wise tij vanishes. We denote the number of sites as N and the 
number of electrons as N e . 

The ground state of the Hubbard model is divided into two 
regions, that is, a weak correlation region and a strong corre¬ 
lation region. When the Coulomb interaction U is as large as 
or greater than the bandwidth 8 t, the ground state should be 
regarded as the strongly correlated one [39, 40]. We use the 
Gutzwiller ansatz and evaluate the superconducting condensa¬ 
tion energy Econd for d-wave pairing. The wave function is 

- Pn^g^bcs, ( 2 ) 

where Pc is the Gutzwiller operator Pq = [1/1 -(1 -g)n^rij{) 
with the parameter g in the range of 0 < g < 1. The gap func¬ 
tion A in the BCS function j/scs is regarded as a variational 
parameter. P^ e is a projection operator that extracts only the 
state with a fixed total electron number N e . The condensation 
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Figure 1: The superconducting condensation energy per site as a function of 
U in units of t for t = 0 and t = 0.2 on 10 x 10 lattice. The number of electrons 
is N e = 88 for t = 0 and N e = 84 for t = 0.2. 


energy is defined as E com i = E(A = 0) - E(A = A opt ) for the 
optimized gap function A opt . 

We show Econd as a function of U in Fig.l. E conc i increases 
rapidly near U ~ St as going into the strongly correlated region. 
Econd is very small in the weakly correlated region for U < St 
and almost vanishes for t — 0. For such small value of E conc i, 
it is certainly hard to obtain a signal of superconductivity by 
means of numerical methods such as the quantum Monte Carlo 
method. It follows from Fig.l that high temperature supercon¬ 
ductivity is possible only in the strongly correlated region. In 
particular, U/t ~ 10-14 is favorable for superconductivity. 
Similar results have been reported for t/t - 0.1, 0.25 and 0.4 in 
Ref. [26], 

3. Enhancement of Superconductivity with Weak Attrac¬ 
tive Interactions 

Fet us examine when the critical temperature T c is enhanced 
in correlated-electron systems. A layered crystal structure is 
obviously an important factor. The increase of the density of 
states in a layered crystal is certainly important and an inter¬ 
layer interaction plays a significant role in finding new super¬ 
conductors [41-45], Apart from this, if there is an interaction 
that promotes superconductivity cooperating with the Coulomb 
interaction, T c will be increased. The d -wave pairing state in 
the Hubbard model will be more stabilized by weak nearest- 
neighbor attractive interactions. We consider the exchange in¬ 
teraction J and attractive interaction V given as 

H = X tijclcjo- + u X n 'T n 4 + J Yj p ' S J ~ \ Hin i 

ij i < ‘j > ' 

+ V ^j n,n i' ( 3 ) 

<<7> 



Figure 2: The superconducting condensation energy per site as a function of A 
in units of t. From the top, we set J = 0, 0.1 and 0.3 for parameters U = 12 and 
t = 0.2. The number of electron is i N e = 84 on 10 x 10 lattice. 


This is so called the t-U-J-V model. In high-temperature cuprates, 
the exchange coupling J arises as a super-exchange interaction 
between d electrons in copper atoms, which is mediated by oxy¬ 
gen atoms. The last term with negative V is expected to arise 
from the electron-phonon interaction [46]. In the variational 
Monte Carlo calculations, both J and V (< 0) act as an attrac¬ 
tive interaction for nearest-neighbour pairs of electrons. The 
Gutzwiller ansatz is also adopted here. We use the electron- 
hole transformation to fix the electron number in tl/ncs instead 
of using /p [47]. We show the results in Figs.2 and 3. In Fig.2 
we compare E coru i for J/t - 0, 0.1 and 0.2 where U/t = 12 
and t/t = 0.2. Basically the exchange coupling J works effec¬ 
tively in a similar way to U to form the d-wave pairing. We 
obtain a similar result for the interaction V, where the lowering 
of the energy E con j by V is smaller than that by J. It should 
be noted, however, that J alone can not stabilize the d-wave 
pairing, which means that if U is small, the d-wave supercon¬ 
ducting state is not realized even for J > 0. This is shown in 
Fig.3 where E con( i is shown as a function A with J - 0.1 and 
t = 0 for U/t — 8 and 12. E con d is very small for U/t — 8 
and vanishes for small U. This indicates that the strong on-site 
correlation is needed so that the exchange coupling J induces 
superconductivity. 

4. Summary 

We have investigated the two-dimensional t-U-J-V Hubbard 
model. The condensation energy becomes large in the strongly 
correlated region, suggesting a possibility of high-temperature 
superconductivity. We mention here that there is a report that 
suggests a high-temperature superconductivity in the moderate- 
17 region by means of the dynamical mean field theory [25], We 
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Figure 3: The superconducting condensation energy per site as a function of 
A in units of t. The condensation energy as a function of A in units of t for J = 
0.1 and t = 0 where U = 8 and U = 12. 


must further clarify a consistency with the results obtained by 
the quantum Monte Carlo method. 

The exchange interaction J and the nearest-neighbor attrac¬ 
tive interaction V cooperate with the Coulomb interaction U 
to enhance superconductivity. It is important that the exchange 
coupling J alone cannot promote superconductivity if U is small 
in the weakly correlated region. We expect that the strong U 
and weak J and V may realize high critical temperature 
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